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THE APPROXIMATE DETERMINATION OF THE FORM 
OF MACLAURIN'S SPHEROID* 

BY 

GEORGE HOWARD DARWIN 

Preface. 

Spherical harmonics render the approximate determination of the figure of a 
rotating mass of liquid a very simple problem. If p be the density, e the ellip- 
ticity, and w the angular velocity of the spheroid, the solution is 

0)2 _ 8 

2^ "IS"''" 

This result is only correct as far as the first power of the ellipticity, but M. 
PoTNCARfi has recently shown f how harmonic analysis may be so used as to give 
results which shall be correct as far as squares of small quantities ; and I have 
myself used his method for the determination of the stability of the pear-shaped 
figure of equilibrium. \ 

Both these papers involved the use of ellipsoidal harmonic analysis, and it 
would be rather tiresome for a reader to extract the method from the complex 
analysis in which it is embedded. It therefore seems worth while to treat the 
well-worn subject of Maclaurin's spheroid as an example of the method in 
question. It will appear below that it would have been possible to obtain a 
more accurate result than that stated above, even if the rigorous solution of the 
problem had been beyond the powers of the mathematician. 

My own personal reason for undertaking this task was that I desired a sort 
of collateral verification of the very complicated analysis needed in the case of 
my previous investigation. 

§ 1. Method of defining the spheroid. 

Let a sphere 8 be described concentrically with the spheroid, and let it be 
sufficiently large to enclose the whole of the spheroid. I call M the region 

* Presented to the Society December 29, 1902. Received for publicatioQ December 13, 1902. 
t Philosophical Transactions of the Royal Society, vol. 198A (1902), pp. 333-373. 
tib., vol. 200 A, pp. 251-314. 
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between the sphere S and the spheroid ; and suppose the density of the liquid 
in xS^ to be + p, and that in -S to be —p. 

If Si denotes any surface spherical harmonic of colatitude and longitude 
()>, it is usual to define the corresponding deformation of a sphere of radius a by 
the equation r = a(l -\- eSl). But in the present investigation it will be found 
that there is a great saving of labor by defining it by the equation 

r^=a\l + 3eSD. 

The two forms give identical results as far as the first power of the ellipticity 
e, but not so when we are to consider the squares of small quantities. 

In general I define S'. by one of the two alternative forms Pl(ix.)^fls<f>, 
where /* = cos 6 . But in the case of the second zonal harmonic (s=0,i=2) 
it is convenient to write 

The fourth zonal harmonic will occur explicitly below, and in accordance with 
the general definition to be adopted we have 

The angular velocity is to be denoted by (o, and the colatitude or cos~V is 
measured from the axis of rotation. 

We must now assume a general form for the equation to the spheroid, and 
shall subsequently determine the several ellipticities so that the surface may be 
a figure of equilibrium. 

The radius of iS being denoted by a, we may write the equation to the sur- 
face of the spheroid in the form 

r' = a' [1 _ 3c + SeS, + SfS, + 3 E/: /^J] • 

In this expression e is the ellipticity corresponding to the second zonal har- 
monic, and it represents that term which exists alone in the ordinary approxi- 
mate solution. Then I suppose thatyandyi are quantities of the order e^ and 
that there are f^ corresponding to all possible harmonics excepting the second 
and fourth zonal ones. Thus all the f. are of order e^, excepting /^ and /^ 
which are zero. Lastly c is an arbitrary constant, and is only subject to the con- 
dition that it is greater than the greatest positive value of eS^ +fS^ + 11 fl S'^ . 
This condition ensures that S shall envelope the whole spheroid. 

It is now convenient to replace the radius vector r by a new variable t , 

defined by 

ff^ — r' 
(1) ^ = -W- 
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Thus the equation to the spheroid may be written 

r = c-e8,-fS,-T.AS\. 

The problem will be solved by making the energy of the system stationary. 
It will therefore be necessary to determine the energy lost in the concentration 
of the spheroid from a condition of infinite dispersion. This will involve the 
use of the formula for the gravity of 8, and since the whole region R is inside 
S , we only require the formula for internal gravity. 

If we were to continue the developments from this point all the formulae 
would involve the constant c . But since it is merely needed for defining a 
sphere of reference of arbitrary size, it cannot finally appear in the formula for 
the energy. It is useless to encumber the analysis by the introduction of a con- 
stant which must disappear in the end, and it is legitimate and much shorter to 
treat c as zero from the first. It is however easier to maintain a clear concep- 
tion of the processes if we continue to discuss the problem as though c were not 
zero, and as if 8 enveloped the whole spheroid. With this explanation we may 
write the equation to the spheroid in the form 

(2) ^=-68,-/8,-^,1181. 

§ 2. The lost energy of the system. 

If the negative density in It were transported along conical tubes emanating 
from the center of 8, it might be deposited as surface density on /S; I refer 
to such a condensation as — C . I do not, however, suppose the condensation 
actually effected, but I imagine the surface of 8 to be coated with equal and 
opposite condensations + C and — C . 

The system of masses forming the spheroid may then be considered as being 
as follows : 

Density + p throughout 8, say -f 8. 

Negative condensation on 8, say — C . 

Positive condensation on 8 and negative volume density — p throughout S, . 
This last forms a double system of zero mass, say D , and Z) = C — R. 

The lost energy of the system clearly involves the lost energy of each of these 
three with itself, and the mutual lost energy of the three taken two and two 
together. Thus the lost energj' may be written symbolically 

\88+lCC+lDD- 8C+ 8D- CD. 

Since Z) is C — R, the last three terms are equivalent to 

- 8C+{8- C){C-R) 8R+ CR- CC. 
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Thus the gravitational lost energy is 

^SS - SR+CR-ICC+ \DI). 

The lost energy of the system, as rendered statical by the imposition of a 
rotation potential, is clearly \Au>^ ^ where A is the moment of inertia of the 
spheroid about the axis of rotation. 

If A^ denotes the moment of inertia of the sphere S, and A^ the moment of 

inertia of the region It considered as being filled with positive density + /a, we 

clearly have 

A=A -A . 

s r 

Thus if ^ denotes the lost energy of the system as rendered statical by the 
imposition of a rotation potential 

(3) B=ISS- SIi+ OR - ^CC + IDD + lci,\A^- a;). 

§ 3. The energy ^SS- SR+ CR-l CC. 

It is in the first place necessary to obtain certain preliminary analytical and 
numerical results. 

If we write da for dfi d(j) , it is clear that an element dv of volume is given by 

dv =■ a? dr da = ^ dr da , 

where Jtf is the mass of the sphere S. 

When we integrate throughout the region R the limits of t are — eS^ — fS^ 
- Z/i 'S'J to zero. 

I now define certain integrals, viz : 

4>\ = i^f(Sirda, a>l = ^JiSJSlda. 



It is well known that 



X 



'+1 2 (i + s)! 



and since in every case, excepting that of /S^, Si = RKfi) liaS4>, we have 

.s L^ (i+i) ! 

^'~2i+l(i-s)!' 

Since S^= — IPj ( /'^ ) » t^® value of <j)^ is derivable from the same general 
formula. Hence we have 



(4) </'o = 3, <t>,^^%, <^, = 1 



?• 



1903] THE FORM OF MACLAURIN'S SPHEROID 117 

Since S^ involves either cos s(}> or sin s<f> , to! vanishes unless 5 = 0; hence we 
need only consider to^. 

The function ( S^Y may be expanded in terms of zonal harmonics. Assume 
then 

Multiplying both sides by SSfji-rr, and integrating throughout angular space, 
we find 

«,= „,>„ and (S,f=j:^^S,. 
But, by actual substitution. 



(^) 


% 4 




<^2 


4 
~21' 


^^4 


8 
35' 


and all the 


higher co's vanish. 












It will be noticed that ft), = 


<^2 


, hence we have 


«o/</'( 


,= i<^. 


(4) 


4 
'"0=15' 


">2 


.= — - 


16 

3^5.7 ' 


«4 = 


8 
■ 3.5.7* 



Also 



The equation to the spheroid possesses a certain property which permits us to 
effect a great saving in the subsequent work. If k be any arbitrary number it 
is clear that to the order of squares of small quantities we may write the equa- 
tion to the spheroid in the form 

since we have only imported a new term of the order e' . Now when k is zero 
the energy will be found to involve terms in e% e', e*, e^f, /"% (fiY- ^ tJ^^n 
we write e + ke^ for e , as is clearly permissible, the term in e^ will give rise to 
a term in 2^e*. Hence it follows that the term in e* in the energy is really 
indeterminate, and that the retention of it would give rise to a fictitious 
accuracy. It is therefore permissible to omit the term in e*, while retaining 
other terms of the fourth order. Again the moment of inertia will involve 
terms in e, e^, e', ef, and the same argument shows that the retention of the 
term in e* would give rise to a fictitious accuracy.* 

It is now necessary to determine the volume of the region i? ; it is 

a'JJdrda = - a' j'\_eS, +fS^ + Tfl S]] da = . 



* If the coefficient ot the terms in «^ and e were zero, the term in e* would not be fictitious. 
This was the case in discussing the pear shaped figure of equilibrium. 
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From this it follows that the mass of S is equal to that of the spheroid, and 
therefore If is the mass of the spheroid. It is this result which makes the 
choice of t as independent variable so convenient. 

We are now in a position to determine the several contributions to the lost 
energy. 

The lost energy of the sphere, denoted by ^SS, is known to be HP/ 5a. 
This is a constant and may be omitted as being of no further interest. The 
internal potential of S is given by 

F=f7rp(3a^-r2). 

But r'=«'(l — 3t), and r^ — a^{l —2t — t^ ■ ■ ■). Therefore as far as 
squares of small quantities. 

Since the' volume of S is zero the first term of V contributes nothing to the 
lost energy SH , and the second term of V will give the whole. Therefore to 
the fourth order 






- s^paf^''^ ^^^' + ^^^'^^^^ + ^^^fi^^^i +r{sj + 21://: ^.^: 

Thus on rearranging the terms we have 

(5) ISS-SR=^ l-Xe^4>^+Xe'a>^ + K/o,,- i/^</.,- J Z (y:)'</>:]. 

We have next to consider the terms depending on the condensation C . Since 
dv/dr = a^dcr the amount of matter in the region S, if filled with density + p, 
which stands on an element of unit area is 

pa JdT ^-pa{eS,+ fS, + L/^ S', ) . 

This expression gives the surface density of the condensation + C, and it is 
expressed in surface harmonics. 

Now by the usual formula of spherical harmonic analysis the internal poten- 
tial of surface density paeS] is 
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2i + l a'--' '~ (2i + l)a «* '" 

As far as the first power of t, 

r' = a'(l — tV), . 

but in the case i = 2 , as far as squares of small quantities, 

»-2=a2(l_2T-T^). 

Hence it follows that the internal potential, say V^ , of the condensation + C is 
given by 

On multiplying this by ^MdTdaj4tTr and integrating throughout R we shall 
obtain the lost energy CR. 
Thus 

crJ-?-((v.cM. 



9Jf2 

47ra 



= 47rJ J 



9J/' 

47ra 



The term in e* may be now omitted in accordance with the principle explained 
above, and therefore 



3Jf ' I , „ . , .. . _ 1 



'^^ = ^ I i^'*^^ + ^■^"'^^ + S 2ii-i (■^*)''^^ 
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In order to find J CO we have only to deal with surface density. Then the 
value of P^ at the surface is given by (6) with t = ; therefore 



*F = 



2''~ 2a 



T 11 



(8) 



An element of mass of the surface density + C is 

Q TIf 

Multiplying these two together and integrating, we find 
And subtracting this from (7) 

+ ii:2^(/:-)v:]. 

Again adding this to (5) we have 

(^) •- ;_i -1 

It remains to determine the value of the term ^DD, and for this end we must 
investigate the theory of double layers, according to the ingenious method de- 
vised by M. Poincar6. 

§ 4. Double layers.* 

Let a closed surface S be intersected at every point by a member of a family 
of curves, and let a be the angle between the curve and the outward normal at 
any point. At every point of S measure along the curve an infinitesimal arc t , 
and let t be a function of the two coordinates which determine position on S. 
The extremities of these arcs define a second surface S', and every element of 
area da of S has its corresponding element dcr' on <S". Suppose that S is 
coated with surface density S, and that /S' is coated with surface density — S', 
where Bdcr = 8' dcr'. The system SS' may then be called a double layer, and 
its total mass is zero. We are to discuss the potential of such a system. 

Let C/^ and U_ be the external and internal potentials of density S on S, 
and f7„ their common value at a point i-* of S. At I* take a system of rec- 

* This section is contained in my paper on the Stability of the Pear-shaped Figure of Equilibrium, 
Philosophical Transactions of the Royal Society, vol. 200A (1903), pp. 251-314 but 
is reproduced here in order that the present investigation may be complete in itself. 
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tangular axes, n being along tlie outward normal, and s and t mutually at right 
angles in the tangent plane. 
In the neighborhood of I* 

dU. dU+ dU. 

dU_ dU_ dU^ 

In the first of these n is necessarily positive, in the second negative. 
Now dU^jds = dU_/ds = dU/ds ; and the like holds for the differentials 
with respect to t. 

Also by Poisson's equation, 

dU_ _ dU^ _ ^^g 
dn dn 

Let PP' be one of the family of curves whereby the double layer is defined, 
and let P' lie on S', so that PP' is t. By the definition of a the normal 
elevation of S ' above /S' is t cos a . 

Let V, v' be the potentials of the double layer at P and at P' . 

The potential of S' at P' differs infinitely little in magnitude, but is of the 
opposite sign from that of S at P ; it is therefore — U^ . The point P' lies 
on the positive side of /S at a point whose coordinates may be taken to be 

n, = T cos a, s = T sin a, t = . 

Therefore the potential of S at P' is 

dU. . dU 

Cr„ + TCOsa-^^ + Tsma^-. 

Therefore 

dU. . dU 

■y = T COS a — 7 — + T Sin a -^j— ■ 
dn ds 

Again the potential of /S' at 2^ is £7^ , and since P lies on the negative side 
of S' and has coordinates relatively to the n, s, t axes at P' given by 

71 = — T cos a, s = — T sin a, t = 0\ 

since further the density on *S" is negative, we have 



dl7_ . dU 

V = T COS a 

Therefore 



V = T COS a — ^ 1- T sin a —,- 

dn ds 



V — V = T COS 



a —, , — = 47rTO COS a . 

|_ aw dn J 
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The differential with respect to n of the potential of S falls abruptly by 47rS 
as we cross S normally from the negative to the positive side ; and the differ- 
ential of the potential of S' rises abruptly by the same amount as we pass on 
across ,tS" . It follows that dv/dn on the inside of S is continuous with its 
/alue on the outside oi S' . 

The surface S to which this theorem is to be applied is a slightly deformed 
sphere, and the curves are radii drawn from the center of the sphere which is 
deformed. The radii are normal to the sphere, and where they meet ^S" the 
angle a will be proportional to the deformation whereby S is derived from the 
sphere. It follows that cos a will only differ from unity by a term propor- 
tional to the square of the deformation, and as it is only necessary to retain 
terms of the order of the first power of the deformation, we may treat cos a 
as unity. 

We thus have the result 

V — v' ^ 4ttt8 . 

Suppose the curve J'l'' produced both ways, and that Jf^ , M^ are two points 
on it, either both on the same side or on opposite sides of the double layer. 

Let M^M^ be equal to f, let f be measured in the same direction as n, and 
let f be a small quantity whose first power is to be retained in the results. 

Let V(|, «j be the potential of the double layer at M^ and M^ respect- 
ively. 

When f does not cut the layer we have 

dv 

and when it does cut the layer 

In the application which I shall make of this result the surface S' will actu- 
ally be inside S. Then v^ will denote the potential at any point not lying in 
the infinitely small space between 8 and S' , and v^ is the potential at a point 
more towards the inside of the sphere by a distance f ; S is the surface density 
on the external surface 8 and t is measured inwards. If then we still choose 
to measure n outwards, as I shall do, our formula becomes 

«„ — «, — K-,— = 47rT8 or , 
" ^ dn 

according as f does or does not cut the double layer. 

It may be well to remark that v being proportional to tS, ^dv/dn is small 
compared with iirrS . It is also important to notice that the term 477x8 is inde- 
pendent of the form of the surface, and that dv/dn will be the same to the first 
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order of small quantities for a slightly deformed sphere as for the sphere 
itself. 

We have now to apply these results to our problem. 

The position of a point in the region R may be defined by the distance meas- 
ured inwards from the sphere 8 along one of the radii orthogonal to S. The 
surface of the spheroid as defined in this way is given by e, a function of 6 and 
4> • Any point on a radius may then be defined by se , where s is a proper frac- 
tion. If s is the same at every point the surface s is a deformed sphere ; s = 1 
gives the spheroid and s = the sphere S. 

If d(T is an element of area of S , the corresponding element on the surface s 
will be ( 1 — Xes ) da- . The value of X will be determined hereafter, and it is 
only necessary to remark that it is positive because the areas must decrease as 
we travel inwards. 

Let s and s + ds be two adjacent surfaces ; then the mass of negative density 
enclosed between them in the tube of which 

( 1 — Xes ) da- and [ 1 — Xe ( s -f cZs ) ] cZo- 

are the ends is — pe ( 1 — Xes ) da ds . If this element of mass be regarded as 
surface density on s , that surface density is clearly — peds. If the same ele- 
ment of mass were carried along the orthogonal tube and deposited as surface 
density on S that surface density would be — /3e(l — Xes). The sum for aU 
values of s of all such transportals would constitute the condensation — C 
already considered. 

The double system D consists of the volume density — p in H, and the posi- 
tive condensation + C on S, the total mass being zero. 

Let z, a proper fraction, define a surface between the sphere S and the 
spheroid. Consider one of the orthogonal lines, and let F„ be the potential of 
D at the point R where the line leaves S and V^ the potential at the point Q 
where it cuts z . Then I require to find V„— V^. 

Since s denotes a surface intermediate between S and the spheroid, 
dsd(^Vg — V^)/ds is the excess of the potential at P above that at Q of surface 
density — peds on s and surface density -f /3e(l — \es)ds on S. Such a sys- 
tem is a double layer, but there is a finite distance between the two surfaces, 
and the form of d(V^— VJ/ds will clearly be different according as z is 
greater or less than s . 

The arc es may be equally divided by a large number of surfaces, and we may 
take t to define any one of them. Now we may clothe each intermediate sur- 
face t with equal and opposite surface densities ±pe[l — Xe(s — t)'\dt. 

The density 

+ pe [ 1 — Xe ( s — ^ ) ] fZ< on t, 

together with 
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— pe^l — Xe^s — t — dt)'jdt on t + dt, 

constitutes an infinitesimal double layer ; and since the positive density on each 
t surface may be coupled with the negative density on the next interior surface, 
the finite double layer may be built up from a number of infinitesimal double 
layers. Hence dtdtd^{V^ — V^)/ds dt is the excess of the potential at P above 
that at Q of an infinitesimal double layer of thickness edt , and with surface 
density pe [ 1 — Xe ( s — t)^ dt on its exterior surface. 

We may now apply the result v^^— v^ — ^dv/dn = 47rST or , according as f 
does or does not cut the double layer, and it is clear that 



d" 
dsdt 



r^[n-K-^«^] = 4,rpe^[l-^<«-0] or 0, 



according as a is greater or less than t. 

In the next place, we must integrate this from t = s to t = , and the result 
will have two forms. 

First, suppose a > s ; then for all the values of t, zy> t, and the first alter- 
native holds good. Therefore 



(^o-T^.-^«|^)=47rp^^[«-i^^«^] 



Secondly/, suppose a < s ; then from t = s to t = z, z -Ct and the second 
alternative holds, while from i = ato^ = 0,a>^ and the first holds. There- 
fore 

ds ' 



:[v,-V^-ez^^^ = 4^pe^[z-Xe{sz-iz^):\. 



We have now to integrate again from s = 1 to s =. . 

From s = l to s = a,a<s and the second form is applicable ; from s = z to 
s = , a > s and the first form applies. 
Therefore 






4.Trpe^ I [ a — Xe ( sa — Ja^ ) ] (?s -f 47rpe- I [ s — JXes^ ] ds 

Jz Jo 

= 47r/>e={a(l -a) - X€ [ia(l - a^) - ia^(l - a)] 

+ ia^-lXa'} 
= 2Trpe^ { 2a - a^ - Xe(a - a^ + Ja=)} . 

Finally, we have to multiply — J ( F^ — K ) % ^° element of negative mass 
at the point defined by z and integrate throughout R. The physical meaning 
of this integral will be considered subsequently. 
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We have already seen that such an element of mass is given by 

— pdv = — /3€ ( 1 — Xes ) d<y dz 

and the limits of integration are a = 1 to s = . 
Therefore 

^ f{V,- VJpdv = 7r/3= r Ce'{l - \€z){2z - z' - Xe(» - z' + lz^)}dzda 

r r dv 

+ iP I I i^^Ci- — '^ez)—j— dzda-. 



In this expression we neglect terms of the order e° and note that e^z^dV/dn 
is of that order. 
Thus 



W ( F„ - V,)pdv = V r r«'[2s - s^ _ Xe(s + »« _ |s')] dzda- 

c r dv 

+ Jp I \ e^z~^— dzda (2 = 1 to 0), 



the integrals being taken all over the surface of the sphere. 
We must now consider the meaning of the integral 



j{V,-V^)pdv. 



Let ^ be a point on S and Q a point on li on the same orthogonal 
line. 

Let — U he the potential at Q of the density — p throughout R , and — C7^ 
its value at P . 

Let S be the surface density of the positive concentration on S, TF its poten- 
tial at Q, and W^ its value at P. 

The lost energy of the double system consisting of — p throughout i? , and S 
on S is 



J r Updv + J r W, Bda - 1 C U^Bda - ^ C Wpdv . 

This is equal to 

if(U- W)pdv - ijiU, - TFJJ 



)hda 
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Consider the triple integral 

Here dv = e^l — Xes )dcrds, U^ — TF„ is not a function of 5 , and the limits of 
s are 1 to zero. Therefore 



But 



t 

I e ( 1 — Xes ) pds 
is equal to h the surface density of concentration. Therefore 

//[C^„ - W,-]Ma = j jj{U, - W,)pdv. 

We may now revert to the Gaussian notation with single integral sign, and 
we see that the lost energy of the system is 



j\_{W,-U,)-{W-U)-\pdv. 



But TF— C^ is the potential of the double system at Q, and is therefore V^% 
and TFp — U^ is the potential of the double system at P , and is therefore FJ, . 
Accordingly the lost energy 

\DD^\j{V,-V, )pdv = l-rrp' J( e' - Xe* )d<T + Ip Je^ ^^ da . 

§ 5. The energy \DD. 

The element of surface of the sphere was written da- in § 4, but in order to 
accord with the notation used elsewhere we must now write it <j^ da- . 
The first term in \DD was 

l-wp^ Ui -\i^)da , 
and when the notation for the element of surface is changed we may write it 

i*/(.-x...... 
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In this expression e is the length measured along a radius from the sphere S 
to the spheroid. We have denoted the outward normal by n , and therefore to 
the second order of small quantities, 

— dn= — dr= a{l + 2t ) Jt . 

The distance measured inward from the sphere to the point defined by t in 
the region i? is, to the first order of small quantities, — n= ar. Again 

e=— I c?w= a I (1 + 2t)Jt 
(10) 

Since — n is what was denoted es in the general investigation of § 4, we have 

(?» = — a^(l -{- \n)d7i dcr 

= a3(l-XaT)(l + 2T)dTda- 

= a'[l + {2-Xa)T'\dTda-. 

But since dv = a^drdcr, we have X = 2/a. 
Therefore 

^=-a^le\SJ+'6eyiSJS,+ STeTc{SJS'i-^e\S,y], 

Xe* = 2a'e\S,y. 

Whence 

This must be multiplied by ^Mpja and integrated throughout angular space. 
Then since, as before, we may omit the term in e'*, this contribution to the 
energy becomes 

(11) ~[-K»,-|eya,J. 

The second term 

dV 



i"/- 



, d(T 
dn 



in \DD remains for consideration. It will clearly be a term in e* and as such 
might be omitted, but it is of some interest to see how it may be computed, and 
I therefore proceed. 

In order to evaluate d Vjdn it suffices to imagine the volume density — p in 
the region It concentrated on a surface bisecting the space between 8 and the 
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spheroid. We may then treat the system D as an infinitesimal double layer of 
thickness \e and with density + C or — pa{eS^ +f>Si + ^.f\^\ ) on its outer 
surface. In the present instance it suffices to consider only the leading terms 
in the density and thickness. Hence by (10) the product denoted t8 in the 
general investigation becomes 

We thus have tS expanded in surface harmonics. 
Now consider two functions 

F^ = ^ A.—^:^ S'^, for space external to S, 

Vi= —^ -A.^ j^^ — S\, for internal space. 

They are solid harmonics and as such satisfy Laplace's equation throughout 
space. Hence they are the external and internal potentials of a distribution of 
matter on S, but since they are not continuous, while their differentials are 
continuous, that matter constitutes a double layer. 
At the surface r ^ a, 

But this must be equal to 477x8. Hence 

2,rpaV [#,^„ + ^-^ ^^ + ^1 ^^] = rio s^ + ^1 ^^ + ^1 S,. 
Therefore 

A^ = I Mae' • i./., ; A, = %Mae- ■ 1 ^-- ; A, = ^Mae" ■ | ^* . 
Now 

dn ~ dr ~ dr ~ ' 



Then since 



,eV = ,«V(^)^ = ^e^[.<^, + ^^ + J^], 
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we have (on writing ahla for the da- of the general investigation of § 4) 

(12) 

- ^ a U^ <^. ^"'~ </>. J' 

This term involves e*, and may be omitted, but, as stated above, it seemed 
worth while to show how it may be computed. 

Then adding (11) to (9) and omitting (12), we have for the whole gravitational 
lost energy 

If 2 r ; _ 1 n 

(1^) ^ L " ^'^'^■' "^ ^'^''"^ ^ ^*''-^"^ " ^-^'"^^ " ^ 2IT1 ^^'^ )'^^ J ■ 

§ 6. Moment of inertia. 

Since ft)^ is of order e , the moment of inertia must be determined to the cubes 
of small quantities. 

We have, to the square of t, 

x^ + 2/^ = r=sin^^ = ^^(|+^,)(l-2T-T'). 

The region R is to be considered as filled with density + p , and the element 
of mass is ZMdTdcrjA^ tt. 

Hence the moment of inertia of the region li is 

' \4:T7 J pa 

+ 2efS, 8,+ 2^ ef\ S, S^ - h\SJ] drd^r . 
The term in e' may be omitted for the reason assigned above, and we have 

+ ^i:efl{8JSl+ie^'{SJ]da 

The moment of inertia of S is | J/a^ . Therefore 

^'~ iirpa^^'' 

Trans. Am. Math. Soc. 9 



Jf{i+S,)(l-2T-T^)drda, 
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Then S M^ 

C=e.-C. = 4— [f + (e + fe^)<^, + e^a,,+ 2e/a,J . 

Whence 

(14) i W = ^' £l [3 + ( 3e + e^)^, + fe^o,, + 3e>J . 

§ 7. Solution of the problem. 
from (13) and (14) we have 

^^ -^ = - ¥<!>. + 1:*5«'«. + \¥'f<^. - ir<t>. - Z 2rri (/:)'<^: 

+ ^ [f + (l« + «')<^2 + K«2 + 3e>J. 
The conditions for a figure of equilibrium are 

^_n ^-0 -^-n 

with 0)^ constant. 

The last of these gives at once /J = 0, and the two others give 

- h<f>. + ie'<o, + ffe>, + ~[(| + 2e)4>, + 3ea,, + 3>,] = 0, 
(15) 

From the former of these as a first approximation 

_= _4_e 

47r/3 ^5 ■ 

The second equation then gives 

Therefore 

(16) f=¥^ = ieK 

We see that y is of order e^ , as was assumed to be the case. Now it is of 
no use to retain terms in e^ in the first of (15), because we have neglected the 
term in e* in B. Thus the first of (15) reduces to 

- ie4>, + le^o,, + £^ [(I + 2e)<!>, + 3ea,J = 0. 
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Whence w^ 



[l^|...eg] = A.(l-.e»;); 



or 



£7 = l*5«(l-|«-4ej) = ^V[l-i«+H«]. 
(17) 

= tV(1-|-0- 

It follows from (16) and (17) that the equation to the surface of equilibrium is 

(18) 

and 2^ = ^56(1-^6). 

It remains to verify that the solution (18) is correct. 

The equation to an ellipsoid of revolution, whose equatorial and polar radii 
are a^ and a^{l — e^) , is 



a\ 



; + sin= 



If we determine r' by developing this expression as far as e\ , it appears that 
the equation to the ellipsoid may be written in the form 

Since the volume of this ellipsoid is |7ra'(l — e^) and that of our spheroid of 
equilibrium was ^ira? it follows that 

a* = Oj ( 1 — e, ) . 
If then we write 

« = «i + Tf «? ' e^ = cj , 

the equation to the ellipsoid of revolution becomes 

and this is the form determined above in (18). 

If the eccentricity of the ellipsoid be denoted by sin 7, we have 

003^7 = (1 - eJ2 = ( 1 - e + ^^fe^)^ = 1 - 2e + Jfe' • 
Therefore 

cos 7 = 1 — e + -j^e^ , 

sin- 7= 2e(l — |e) ; 
whence 
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C0S7 sin^7 = 2e(l — ^-e), 

cos 7 sin* 7 ^ 4e^ . 

Now it is known that the rigorous solution for the angular velocity of Mac- 
laurin's ellipsoid may be written in the form 

ffl^ _ ^ (2w- l)!sin2''7 

^ _ cos 7 2. (2m + 1)(2«, + 3) [{n-iy.f'A^"-^ ' 

Taking the first two terms 

4^3 = i\ sin=7 cos 7 + -^^ sm*7 cos 7 

This agrees with the second of (18), and the solution is found to be correct 
as far as squares of small quantities. 

The approximate solution found above is insufficient to enable us to discuss 
the stability of the Maclaurin figure, but it may be well to indicate how a more 
accurate approximate solution would give the required residt. 

Let us suppose that V is the gravitational lost energy corresponding to the 
equilibrium eUipticities e, f, /I and the angular velocity w . 

Then we are to regard V and the moment of inertia C as functions of 
e , /', f\ , so that e , /' , f\ are the solutions of the equations 

dV dC dV dC dV dC 

Now suppose that the eUipticities corresponding to any neighboring form are 
e + Be , f + Bf, fl + Bf. , and that V+ BV, C + BC are the corresponding 
values of V and 0. Conceive also that these variations from the equilibrium 
figure are made subject to constancy of the angular momentum, then it is clear 
that IT, the sum of the potential and kinetic energies, is given by 

U=-iV+BV) + i-^^^^. 

If we omit constant terms and only retain squares of small quantities we have 

U= -BV- i»-S(7 + i»2 ^^ . 

For the sake of brevity I will only retain the two eUipticities e , f, since this 
will suffice to indicate the general law. Now if V and C be expanded by 
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Taylor's theorem in powers of Be, Sf, and if we bear in mind the conditions 
(19) for the figure of equilibrium, it is easy to how that 

Hence 

This is a quadratic function of 8e , Bf and the vanishing of the Hessian would 
give the condition for the change from stability to instability. However, we 
need just one term more than that found above to obtain even a first approxi- 
mation to the limit of the stability of the Maclaurin spheroid. It is accord- 
ingly useless to pursue the topic further. 

Cambbidge, England. 



